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THE PREDICTION OF 

NORMAL FORCE AND ROLLING MOMENT COEFFICIENTS 
FOR A SPINNING WING 


B. W. McCormick 


SUMMARY 


The integrated normal force and rolling moment coefficients 
for a spinning wing are predicted by means of a single strip 
analysis. Non-linear airfoil section data for angles of attack 
from 0® to 180“ are used in a small computer code to nxmerically 
integrate the section normal force coefficients along the span 
as a function of the local velocity and angle-of-attack resulting 
from the combined spinning and descending motion. A correction 
is developed to account for the radial pressure gradient in the 
separated, rotating flow region above the wing. This correction 
is found to be necessary in order to obtain agreement, both 
in form and magnitude, with rotary balance test data. 



THE PREDICTION OF 

NORMAL FORCE AND ROLLING MOMENT COEFFICIENTS 
FOR A SPINNING WING 


B. W. McCormick 
INTRODUCTION 

The purpose of this study is to predict the normal force 
and rolling moment coefficients for a spinning wing. These 
predictions are then compared with experimental measurements 
obtained with a rotary balance. The predictions are made by 
calculating the local angles of attack and local dynamic 
pressures. These determine the wing section normal forces which 
are then integrated numerically to obtain the normal forces and 
rolling moment coefficients for the entire wing. 

It will be shown that a simple strip analysis, as just 
described, underestimates the magnitude of the wing’s normal 
force coefficient and also fails, qualitatively, to predict the 
effect of spin rate on Cj^. The strip analysis, of course, 
neglects 3-D effects caused by induced velocities from the trailing 
vortex system. However, such effects will generally reduce the 
Cjg prediction and are, thus, not the explanation for the dis- 
crepancy between strip analysis and experiment. After some 
study, it was concluded that the stalled flow above a spinning 
wing results in centrifugal pimping and, hence, a spanwise 
pressure gradient, which produces a significant increment in Cj^. 

A relatively simple expression for this increment as a function 
of spin rate is developed which agrees reasonably well with 
experiment. 



NOMENCLATURE 



wing span, m 
chord, m 

rolling moment coefficient 
maximtim wing lift coefficient 

normal force coefficient 

parameter in expression for (eq. 4c) 

parameter in expression for (eq. 4b) 

slope of Cjj vs ot curve for small a (eq. 4a) 

local static pressure. Pa. 

static pressure in separated flow region 

local dynamic pressure. Pa. 
vertical velocity of descent, m/s 
dimensionless spanwise location, y/ (b/2) 
value of X, outboard of which, wing is unstalled 


y = spanwise distance from wing centerline, m 
a = section angle of attack, degs 

= section angle of attack for left wing, degs. 

= section angle of attack for right wing, degs. 

- parameter defined in figure 3 


a = parameter defined in figure 3 
max 


= denotes increment in 

ACjj = denotes increment in 

AN = denotes increment in normal force 


E = exponent for Cjj curve fitting (eq. 4b) 



>3 

6 “ angle between chord line and vertical (fig. 1) 

(|) = angle of chord line from vertical, (eq. 3) 

(j = dimensionless spin rate, nb/2V 

n = spin rate, radians/ sec. 

ANALYSIS 

In a steady spin, the right and left wings of an airplane 
are generally operating at different angles of attack. For a 
spin to the right, the resultant flow over the right wing may 
actually be coming from the direction of the trailing edge, 
which in helicopter rotor terminology is designated as reversed 
flow. 

Consider the sketch given in figure 1 which depicts right 
and left wing sections in a spin to the right. The figure on 
the left is a view looking inboard along the left wing and the 
figure on the right is a similar view for the right wing. The 
angle of attack of the left wing section, cl , will be given by, 
ul = 0 - <l> (1) 

and for the right wing, 

“r. = tt - 0 - (j) (2) 

with measured from the trailing edge . 

0 is the angle of the chord line from the vertical while 
the angle (|) is given by, 

V 

n being the angular velocity about the spin axis, y the 
spanwise distance from the centerline and V the descent velocity. 
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If we let, 

X = y 
TbT^ 

- 

n 2V 

then 4> becomes, 

^ = tan”^ Cx (3) 

An appreciation for “R can be gained from figure 2 

which was calculated on the basis of the preceding equations. 

In this figure, contours of constant ul and ur are presented as 

a function of 0 and ilJx . Spotted on the figure are steep and 

flat spin (or moderately flat) test points for the spin tunnel 

model and full-scale version of the typical, low-wing general 

aviation airplane. These points are for the tips with x= 1.0. 

Thus the angles of attack inboard of the tips will be greater 

than those read from the chart. For the model steep spin, 

equals approximately 21° and ur , relative to the leading edge 

is approximately 50°. The corresponding angles for the full- 

scale airplane are only 3 to 5° higher. For the flat spin mode 

however, the angle of attack for the right wing is significantly 

different between the model and full-scale tests. For the model 

tests the flow on the right wing tip is reversed at a value of 

about 60°. On the full-scale airplane, the angle becomes 78° 

relative to the leading edge. 

Since static tests show the model wing to stall at an a 

of about 10° (Ct. 1.0), and even allowing for an increase in 

"max 

CLmax with Reynolds number, it would appear that both wings 
are completely stalled at the test point conditions. 
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DEVELOPMENT OF PREDICTION METHOD 

A small computer program has been developed which performs 
a strip analysis and numerically integrates along the span to 
predict a spinning wing normal force coefficient, and 

rolling moment coefficient, . The form of the static normal 

force coefficient vs. a curve used in the program is shown 

1 

in figure 3, i 


Cn “ CjicjO 



“ ^ “max 

(4a) 

Cn = ‘^Nmax 

(sina) 

'c 

“ > “2 

(4b) 

Cn = - 

na - 

®max“| 

<=>ll - CNj) Vx - “ ^ “2 

(4c) 

- 

“max-1 



This form for the curve was chosen on the basis of 0012 

airfoil and flat plate data. In order to account approximately 

for 3-D effect, the 0012 data was modified as shown in figure 4. 

The slope was corrected for aspect ratio in the usual manner 

while C« was reduced on the basis of the variation with aspect 
•^max 

ratio of the drag of flat plates normal to the flow. The "section" 
curve which was finally used is the lower one in figure 4 
labeled "3-D wlhgi". 

It is interesting to note that the curve is nearly symmetrical 
about an a of 90°. Thus an angle of attack relative to the 
trailing edge will produce the same as that for the same 
angle relative to the leading edge. Therefore, for 6 values 
near 90°, where the angle of attack and resultant velocities 
are nearly equal on the left and right wings, but in opposite 
directions, one would expect the resultant force on the two 
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sides to be equal. Thus the rolling moment should approach 
zero, regardless of S , as 8 approaches 90°. 

In terms of the section values, the wing and can 
be determined as follows. 


N = / 1/2 p[v^ + (ny)^]cNCydy 

./j. 


^b/2 
b/2 


1/2 pCv-" +(ny)'‘]Cu cy dy 

li 



/2 

1 / 2 p[v^ + (ny)^]CNj^Cy dy 


In dimensionless form, for a rectangular wing, these become, 


N = y* [l + 

rir 


= 1/4 


dx 




1 + ({5?i)^|(C 


Nl 


- c 


Nr) 


3? dx 


(5) 

(6) 


The Cjj notation is somewhat confusing. Under the integral 
sign it refers to section data while otherwise, it denotes the 
normal force coefficient of the entire spinning wing. 


INITIAL RESULTS 

Equations (5) and (6) were evaluated numerically using 
figure 4 and eqs. (1) through (4) . The results of those cal- 
culations are presented in figures 5 and 6 together with rotary 
balance data. This data was taken from the reference by sub- 
tracting fuselage-only data from data for the wing and fuselage 


combination. 
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The lack of agreement between predictions and experiment in 

these figures, particularly for is obvious. Of most concern 

is the disagreement in at the higher 9 and <5 values where 

flat or moderately flat spins occur. Varying parameters such 

as C„ , C„ , 01 , oroto within limits did not materially affect 

N 2 max z 

the C„ predictions. Increasing C merely shifts the curves 

max 

up in figure 5 without changing the rate of increase of with 
3 to any appreciable extent. In view of this rate being signi- 
ficantly higher for the data than for the prediction, it was 
felt that something basic was lacking in the analysis. The 
usual induced effects were not the answer since downwash 
corrections would tend to reduce the calculated C„ values even 
further. Also, compared to a helicopter rotor, the "disc" 
loading of a spinning airplane is low and its rate of descent 
high so that its downwash velocities should be low. 

After further consideration of this problem, it was realized 
that the fluid mechanics of the completely-stalled, spinning 
wing was not being adequately modeled. A correction to 
which was developed to account for the rotating separated flow 
will follow. 

RADIAL PRESSURE GRADIENT CORRECTION TO C„ 

' " 

Immediately above a spinning wing, and for some distance 
beyond, the flow is separated. This means that the fluid in 
the separated region is moving with the wing in the form of a 
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solid body rotation. Hence, a radial pressure gradient must 
exist through this fluid given by. 


Integrating, this becomes, 

p + const. (8) 


At r=b/2, p must equal the pressure p^ along the separated 
streamline at the tip. Thus, using the point to determine the 
constant in equation (8) results in, 

P - Pc = “ ^ (1 - x:2) ( 9 ) 

Without the rotation, p^ would be the constant pressure along 

the back surface of the wing. Hence the pressure decrement, 

p-p , results in a increment to calculated purely on the 

basis of static section normal force coefficients. The 

increment in the wing's normal force can be written as, 

2 

^(Pc - P)dy 

or in coefficient form. 




(10) 


This correction of course only applies to separated flows 
and is independent of 0 provided 0 is sufficiently high for a 
given u , to assure that the wing is stalled. 

Referring to figure 2, only the right wing will be completely 
stalled for a 0 of 30°. Thud, for this case, only half of 
equation (10) was used in calculating 


A. 
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Predictions, with and without the correction for the radial 
pressure gradient are presented in figure 7. The correction is 
seen to have a significant effect on the shape of the curves, 
appreciably increasing the values at the higher 3 values. 

Figure 8 compares the experimental values of with 
predicted values including the correction for the radial 
pressure gradient. Here, the agreement is seen to be con- 
siderably improved over that of figure 5. On the basis of 
this figure and the earlier physical reason, it would appear 
that any calculation of should include the increment given 
by equation (10) . 

The radial pressure gradient should not effect the rolling 
moment, provided both wings are stalled. However, should one 
wing not be stalled, such as the left wing in a right spin, or 
vice versa, then a decrement for a right spin or an increment 
for a left spin should be added to equation (6) . Integrating 
the increment contributed by P-P^ (eg. 9) over one wing and 
non-dimensionalizing results in, 

AC£ = ±a^/16 (11) 

For a right spin this correction is negative. For the 
values of 0 shown here, the left wing is installed over at 
least the outer 70% of its span for a 0 of 30°. At 50° there 
may be a slight correction. Thus the correction given by (11) 
was applied to the calculated values only for 0 =30° and is 
included in figure 6. 
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As a refinement; the corrections to and can be derived 
for the more general case where the advancing wing (left wing 
in a right spin) is unstalled outboard of a station Xs- This 
is done by integrating eg. (8) only out to Xg and setting p 
equal to p^ at that location. In this case equations (10) and 
(ilj become, 

ACn = ^ (1 + Xg^) (12) 

^ (1 - Xg2)2 (13) 

The program appended to this report includes these corrections . 
The program checks to see if any portion of the wing is 
unstalled and sets Xg accordingly. 

CONCLUSIONS 

An approximate method has been developed for predicting normal 
force arid rolling moment coefficients for a spinning wing. Fair 
agreement is obtained with experimentally-measured coefficients 
but only if a correction is included to account for the radial 
pressure gradient produced above the wing by the rotating 
separated flow. 
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LEFT SIDE VIEW 




RIGHT SIDE VIEW 


FIGURE 1 

Velocities Influencing Left and Right 
Sections of a Wing in a Right Spin 
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FIGURE 2 

Local Angles of Attack on Left and Right Wings 
for a Right Spin. (Test points for Wing tips x = 1) 






FIGURE 3 

Model of Curve 
for Airfoil Section 





Airfoil and Wing Normal 
Force Coefficients 

(static) ■ 



FIGURE 5 

Calculated and Measured Normal 
Force Coefficients for a Spinning Wing 




FIGURE 6 

Predicted and Measured Rolling Moment 
Coefficients for a Spinning Wing 



Effect of Radial Pressure Gradient on 
Predicted Normal Force Coefficients for 
a Spinning Wing 
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APPENDIX 


PROGRAM CNCL3D(INPUTjQUTP.UT>TAP13-IN«>UT^.TAP-E6-DliimTl 

THIS_ PROGRAM _PEPFORHS„A_STRIP _ANALYS IL.TQ, X ALCilLAI£._NDRMAL JDRCE 

AND POLL ING MOM ENT C DEFFICIPNTS ON A . WING. 

CNW«W1NG NORMAL FORCE COEFFICIENT 

CLW-WING ROLLI NG MOME NT COEFFICIENT 

CNWC-WING NORMAL FORCE COEFFI CIENT CORR ECT ED FDR RADIAL GRA DIENT 

C L WC - WI NG ROLL ING_MQM£ljTL_CQ£££JLCIENT_XGRR£C TED F QR.-RAJ>LAL_GRiLD-I£]iL 
LIFT CURVE INCREASES .LIN£ARLX_U.P_J0_ANGL£_ALL£_lEQR_L£ADUi6_ED.G.EJ_ 

OF ALTE (FOP TP AI LING f DGEl.. .JT_TiiEN_JlECR£AS-E.S-J.XN£ARlY T O -AN 

ANGLE AL2 WHICH lS_ON-JH£..JCOil£i£.TELY . SJALLEiL-CmVE^ ALJQNG_LH1S 

PORTION OF THE LI FT. CUR VE CN IS AP£RDXIHATEI)__BY_THE JIELATXDNGHIP. 

USED IN THE PROGRAM. TWO ,L INES _BEf ORE. STATE MENJ-Z^ 

FOF A GIVEN.INPUT.J .JCALGUIATIDNS BEGIN , AT.. AhLANGLE . D£_AIIACK_Df- . 

3C DEGFEES AND GO UP_ JD_90. DEGREES, .AT. EACH.ANGLE. C.0MPLU.AII0M.. .. 
ARE P E R F DRMED_I]V£P_A_RANGE_-D F. DIMEN5 1 ON LE S S_S P. IN_,RA T£S. £ROM,-D. IQ. I 


EXPRESSEIL.A.S 
CURV£_5LDEE_ 
CLA-0.072 


_ERAXIUlN_ll£_JQi£.. S EM I -SPAN 
AllE-i e.O ALTE»1A.0 A L 2«30. 


AL L AN GLES A ND LIFT 

DELX..05 CNHAX»1 .2 


JJdE SE ARE. 2-D VALU ES AN Q_GJiOULD BE USE D FOR TH E 


_EX£GJi«t622_ 

STRIP ANALYSIS.^_ 

..Ihi£IA_-LS_RE L ATI.y£--TJ3JWIUG_-Z£RQ LIFT-HNL 


.RE AJ31 5j>.IOO ) C L A^ AL L E > ALTE# AL2^ DEL X>CNMAX >£XPCH 

CHA Y»- 1. 

A L 2iA L 27.5.7 .3 ^ 


ALLE?ALLE_/57,3_ 

ALTE ^A L I £_/ 5 7 .3 

XLA-CLA^7.Z 

CN1LE-_ALLE*C_LA 

C N 1 T E -A L T£ *CLA 

CM?«CNMAX*( SIN(A L2) )**EXPCN 


CNW-G. 

CLW-.0, _ 

□ MEGA;^, 

THETA-3_0_, 

, THETJi«THFTA/57.A_ 



1 CONTINUE 

PMI-ATAN (0MEGA*ABS(X)) 

A L D H A L * f H E_TJ.-P_H I 

A L P H A R « 3 . 1 A 1 19- T H E ILA^. I_ 
VFSOD-l,+(OMEGA+X)*»2 


IF (X, GT.C.JGO. .TO. 2 

ALpwA*ALrHAL 

I F ( ALPWA.LT. ALLE) CN-CLA + ALPHA . _ 

IF ( ALPHA. LT.ALLE)XS»X 

IF (ALPHA. LT,ALL£)CHAY-1, . 

IF ( ALPHA. LT. ALLE) GO . TG 3 . . 

IF (ALPHA.LT.AL2)CN-CNlL£-{ ALPHA-ALLE )7 { AL2-ALL£)+lCNlX£=rH2J 

IF (ALPHA.LT. AL2JGD. TO. 3. _ 

CN-CN.MAX* (SIN(ALPHA) )**EXPCN. .. 

GO TO 3 
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2 CONTINUE . 

ALPHA-ALPHAR . _ 

IF (ALPHA. LT.ALTE)CN»CLA.*ALP_HA 

IF (ALPHA, LT,ALTF)60_TD_3 

IF(ALPHA.LT.AL2) CJN - CKl IE.- LA Lf KA-AL If LZ-LAL 2r-.A Lifi ♦.( C HI J_£Td:N2J 

I F ( A L P H A . L T .-AL21.Q.D-TQ 3 

CK.CNMAX* (SIN( ALPHA ) )**EXPCN 

3 CONTINUf 

D C N W - , !»*VR S QD*CN 

D C L W --DC NW* tJh.* 

. . I F ( X , f 0 . r 1 • )AD„ J 0_& 

C H W - C N W ± ( D C KW± D CUHl) /2a*DEQ 

CLW«CLW-KDCLW4DCLW1 )/2 .*Df LX 

DCNWl-DCNW _ _ 

.. . D C..L W 1 • D C L W j 

X«X + D_ELX_ 

IF ( X,GT,1, )*GD_ TD_5 ' 

go' to I 

A DCNWl- DCNW 

DCLWl -DCLW 

X?.X + D5LX 

GC TO. l 

E CONTINUE 

THETA»THETA*57.3 

IF ( C HAY , E Q.-i •_) X S-1. 0 

CNWC«CNW-fOMEGA»*2/3.*(1. 4XS**3) 

CL WC-CLW-0HEGA**2 /16 »* C l.-XS*»2 ) **2 

W P I T E ( 6 » i 0 1 ) ON F G A ,T HETAy C NW j C L W> CNWC» CL WC 

OMEG_A»gME_GAjt_.2_ 

_ CNW- 0. ' 

.7_rcL.w^d,__ I 

X--1. _ _ __ 

.. THEJA-T.HETAiL5Ja3 _1 

IF (ONEGA. GT. l. )GO TO fe 

GO ..TO ,1 

6 CONTINUE 

C.NEGA-D. 

THETA - THE TA+ 10 . /57.3 

IF(THETA.GT.90./.57. 3)GQ TO 7 

GC TO 1 

7 CONTINUE . ' ~ 

CLA.CLA/57.3 ~ ~ ^ 

ALLE-ALLE.+ 57.3 JIZIZZZ ^ 

ALTt« ALTE+57,3 _ 

AL2-AL2»57,3 _1 ' J 

. WPITE(6,102)CLA,A.LLE»ALTEyAL2>DELX*CN_MAX,EXPCN 

STOP , _ ■ 

IOC .FCIRNAT(7F10,A) . 3 _____ ~ 

101 FCRNAT( 2X., 6HCMEGA-y ^^5, 2y AXj|_6HTHETA«> fb,Zj AX^ 3HCN- / FT.S^jAXji 3JiC 1-y F.7 

1.3yAX,£HCNWC-yF_7.3yAX,5HCLWC.-.>Ffc.A) 

1 0 2 F 0 P N A T ( 2 X y A HX L A^ » F 7 . A » A X , 5HA LL E-> F 7. 2 y A y 5H A LTr«yF7,2>AXyAHAL2*yF7 
1. 2y AXy !rHrEL,X-y_FJ. Ay AXy6HC NMAXj»y F7,2y AXy 6HEXPCN»y F7.3) 

END ' ■ ■ ■ 
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